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Abstract. Recent experimental results on diffractive high Et jet production and diffrac- 
tive open charm production, obtained by the HI and ZEUS experiments at HERA, are 
presented. Various theoretical models for diffractive hard scattering are discussed and 
confronted with the experimental data. 
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1 Introduction 

It is well known that soft hadron-hadron scat- 
tering can be described by Regge phenomenol- 
ogy, which models the total, elastic and single 
and double diffractive dissociation cross sections 
in terms of reggeon and pomeron exchange . At 
high energies, it is the pomeron exchange mech- 
anism that dominates the total cross section 0]. 
However, until recently, our understanding of the 
pomeron in terms of Quantum Chromodynamics 
(QCD) Pi remained fragmentary, at best. 

The observation of events with a large rapid- 
ity gap in the hadronic final state at HERA 
which are attributed to diffractive dissociation of 
(virtual) photons, led to a renewed interest in the 
study of the underlying dynamics of diffraction. 
In deep inelastic ep scattering, the long hadronic 
lifetime of the photon at small Bjorken-a: and its 
variable virtuality allow to make the link with 
diffractive dissociation in hadron-hadron scatter- 
ing, while the presence of one or more hard scales 
enables perturbative calculations in QCD. 

In the framework of QCD, it is tempting to at- 
tribute a partonic structure to the pomeron 
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an approach which proves to be very successful 
in describing various aspects of diffractive virtual 
photon dissociation. QCD calculations indicate, 
however, that the pomeron does not exhibit a uni- 
versal behaviour over the full kinematic range. In- 
deed, experimental observations of a transition do 
exist when going from the soft to the hard scat- 
tering regime, changing the effective intercept of 
the pomeron trajectory from 1.08 in the soft limit 
to 1.2 in harder interactions. 

Early experimental measurements at HERA of 
the inclusive cross section [^|-|^ led to the conclu- 
sion that a partonic pomeron would be dominated 
by hard gluons. This has been qualitatively con- 
firmed by several studies of the properties of the 
hadronic final state of diffractive virtual photon 
dissociation 0, which are all in agreement with 
the fragmentation of a colour octet-octet partonic 
final state, leading to higher particle multiplici- 
ties, lower thrust values and a a final state system 
that is not strictly aligned with the 7* axis in the 
rest frame of the photon dissociation system. 

This article presents some recent experimen- 
tal measurements of more exclusive final statesQ. 
The production of open charm or of jets with 
large transverse energy is particularly sensi- 
tive to the gluon content of the pomeron. More- 



^ Previous results on jet production in diffractive deep 
inelastic scattering can be found in [p^ . 



2 



Pierre Van Mechelen: Final States in Diffraction at HERA 



over, the high Et and large charm quark mass 
provide an additional hard scale for perturba- 
tive QCD calculations. Following this introduc- 
tion, Sec. 1^ provides a short reminder of the kine- 
matics of diffractive deep inelastic scattering. Sec- 
tion 1^ discusses the theoretical models that are 
used for comparison with the experimental re- 
sults. Sec. H presents the preliminary results, ob- 
tained by the HI and ZEUS collaborations, on 
high Et dijet production and open charm pro- 
duction and, to conclude, a summary is given in 
Sec. |. 

2 Kinematics of diffractive ep 
scattering 

The generic diffractive process ep eXY is illus- 
trated in Fig. ^ A photon with virtuality = 
—q^ interacts with the proton to produce the two 
final state systems X and Y , which are sepa- 
rated by a large rapidity gap if their invariant 
masses are small. Cross sections can be expressed 
in terms of and the invariant masses of the 7*p 
system {W), and the X and Y systems (Mx and 
My, resp.). In addition, the following variables 
are usually introduced: 



Xjp 



q- {p-Py] 
q-p 



X 



2q-{p-pY) Q^ + Ml 



(1) 
(2) 



In the proton infinite momentum frame, Xjp is the 
fraction of the beam proton momentum trans- 
ferred to X, while f3 represents the fraction of 
the exchanged momentum carried by the struck 
quark. The squared four-momentum transferred 
at the proton vertex (t) has been neglected in the 
above equations. 

3 Theoretical models 

Theoretical models for diffractive deep inelastic 
scattering range from purely phenomenological 
parameterizations with a large number of degrees 
of freedom to full perturbative QCD calculations. 



X 



p 



largest rapidity 
gap in final state 



Y 




Fig. 1. Illustration of the generic diffractive process 
ep — > eXY. 

The following briefiy discusses those models that 
are used for comparisons to experimental data. 

The partonic pomeron 

When measuring the diffractive cross section in 
deep inelastic scattering at HERA, it was soon 
recognized that the xjp dependence can be mod- 
eled exploiting Regge factorization |1T[ , provided 
that a sub-leading exchange is allowed to con- 
tribute at high values of xjp. In this model, the 
diffractive cross section factorizes into pomeron 
and reggeon flux factors, depending on xjp only, 
and into pomeron and reggeon structure func- 
tions, depending on /? and Q^. The flux factors 
describe the long distance physics at the proton 
vertex and can be parameterized using Regge the- 
ory, while the structure functions are free param- 
eters containing the short distance physics at the 
photon vertex. 

Recently, QCD factorization has been proven 
to hold for the leading twist component of diffrac- 
tive deep inelastic scattering allowing to fur- 
ther interpret the structure functions in terms of 
quarks and gluons by applying a QCD motivated 
model for their /? and dependence Par- 
ton density functions for the pomeron have thus 
been extracted from the measured cross section, 
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assuming the pion structure function for the sub- 
leading reggeon exchange and the next-to-leading 



order DGLAP equations for the evolution in 

The outcome of these combined Regge-QCD 
fits is a value of around 1.2 for the pomeron in- 
tercept aip(O) and a pomeron structure function 
that is dominated by gluons carrying large frac- 
tional momenta W\. 



Soft colour interaction model 

A model that does not use the notion of a pome- 



ron at all, is the soft colour interaction model 
Instead, the standard matrix element and parton 
shower description of deep inelastic scattering is 
used, which, at low Bjorken-x, is dominated by 
boson-gluon fusion. In addition, non-perturbative 
interactions affect the final colour connections be- 
tween partons, while leaving the parton momenta 
unchanged. This may result in an interruption of 
the colour strings between partons, thus creating 
a large rapidity gap in the final state. The proba- 
bility for such a soft colour interaction has to be 
fixed by the experimental data. 

Recently, an improved version of this model 
has become available, which is based on a gener- 
alised area law and which favours soft colour in- 
teractions that make the colour strings between 



partons shorter [16 



Saturation model 

In this and the subsequent models, the virtual 
photon is considered to fluctuate in partonic con- 
figurations, like e.g. qq or qqg. In the proton rest 
frame and at low Bjorken-x, this happens long be- 
fore the actual interaction with the proton. Both 
states, qq and qqg, can be described by dipole 
wave-functions. These wave-functions are respon- 
sible for the basic shape of the /3 dependence of 
the diffractive structure function. Three contribu- 
tions have been introduced UTtI, two of which are 



related to the production of qq pairs by longitu- 
dinally or transversely polarised photons, while a 



third contribution takes the radiation of an addi- 
tional gluon into account for the case of a trans- 
versely polarised photon: 



F^- ~ ^ In 



QQ 



QQ9 



/3(l-/5), 



(3) 
(4) 
(5) 



From these formulae, it is clear that qqg states are 
expected to populate the low (3 (i.e. high mass) 
region, while qq states dominate at intermediate 
and high values of (3. Note that the longitudinal 
qq contribution is of higher twist nature and that 
the longitudinal cross section is expected to dom- 
inate over the transverse cross section at high (3. 

The subsequent models differ in the way they 
treat the actual dipole-proton interaction. 

An attempt to describe the saturation of the 
inclusive cross section at low and low x leads 
to the following parameterization of the dipole- 
proton cross section [0: 



Cr(X,r j = (To 



1 — exp 



with 



1 



GeV \xq J 



mi{x) 



(6) 



(7) 



This dipole cross section exhibits colour trans- 
parency at small dipole radii r (i.e. a ~ for 
r < i?o) and saturation at large r (i.e. a ~ ao for 
r > Rq). Saturation at low x is explicitly built-in 
through the x-dependent radius Rq. 

After a fit to the inclusive structure function, 
the diffractive structure function follows as a pre- 
diction. 

One interesting property of this model is the 
prediction that the ratio of diffractive to inclusive 
cross sections is roughly constant function of 
W , as has been observed experimentally 

Semi-classical model 



In the semi-classical approach [|T^], an approxi- 
mation is used in which the hadronic target is 
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much larger than the transverse separation of the 
dipoles. Assuming further that the colour fields 
in distant parts of the target are uncorrelated, 
the relevant colour field configurations can be av- 
eraged, yielding explicit formulae for both the 
diffractive and inclusive parton distributions at 
some low scale Q^. Diffraction occurs if the emerg- 
ing partonic state is a colour singlet. A conven- 
tional DGLAP analysis is then used to calculate 
the diffractive and inclusive structure functions 
at larger values of Q^. The model has four free 
parameters in total, which are obtained from a 
fit to the experimental data. 

It is interesting to note that also in this model 
the diffractive gluon density turns out to be much 
larger than the diffractive quark density. 

Two-gluon exchange model 

This QCD calculation describes the diffractive 
7*p interaction in terms of partonic fluctuations 
of the photon [qq, qqg, etc.) that scatter elas- 
tically off the proton through the exchange of 
two gluons in a net colour-singlet configuration 
20[| . An essential feature of the model is that 



the cross section is proportional to the square of 
gluon density inside the proton. The calculation 
requires high transverse momentum of all outgo- 
ing partons, thus making it particularly suited 
for diffractive jet production. In order to avoid 
the valence quark region in the proton and to ex- 
clude secondary reggeon exchanges the range of 
vahdity is restricted to xp < 0.01. 

4 Experimental results 

In the following, Monte Carlo models are used 
to compare the measured hadron level cross sec- 
tions to the predictions of the phenomenological 
models and QCD calculations presented in the 
previous section. RAPGAP is used to obtain 
the predictions of the partonic pomeron model 
with different pomeron intercept values and par- 
ton distributions. The programme also contains 
implementations of the saturation and the semi- 
classical models as well as the two-gluon exchange 
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Fig. 2. The differential diffractive dijet cross section 
is displayed as a function of Zp^^ and compared with 
predictions from the partonic ("resolved") pomeron 
model. Different gluon density functions, obtained 
from fits to the inclusive cross section and shown in 
the top figure for an evolution scale which is typical 
for the dijet event sample, are used as input to the 
model. The effect of only including the direct photon 
contribution is also shown. 

4 < Q2 < 80 GeV^ 
0.1 < y < 0.7 



xp < 0.05 
My < 1.6 GeV 
\t\ < 1 GeV^ 
p^'* > 4 GeV 
-3 < r?^'^* < 



Table 1. Kinematic domain of the diffractive dijet 
cross section measurement by HI. 



model. Both versions of the soft colour interaction 
model are implemented in the LEPTO genera- 



tor 22 . 
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Fig. 3. The differential dijet cross section is displayed as a function of xjp (a) and f3 (b) and compared to 
the partonic pomeron model. Predictions are shown for different values of the pomeron intercept ajp(O) and 
also the sub-leading reggeon component is displayed. 



Diffractive high Et jet production 

The production of jets with large transverse en- 
ergy provides an ideal test of the underlying dy- 
namics of diffraction models . The interaction 
mechanism, which is dominated by boson-gluon 
fusion, is highly sensitive to the role of gluons and 
the large transverse energy lends a hard scale for 
perturbative calculations. 

The HI Collaboration recently obtained pre- 
liminary results based on an integrated luminos- 
ity of 17.9 pb~^, resulting in approximately 2500 
events with at least two jets in the final state 
system X. Diffractive events are selected by ex- 
cluding any hadronic activity in the range 3.2 < 
^lab ^ y g jg^g g^j^g identified as collimated 
energy depositions in r/-^ space within a cone of 
radius R = a/ Arj'^ + Acj)^ = 1. The fraction of the 
pomeron momentum transferred to the dijet sys- 
tem, and therefore entering the hard scattering, 
can be calculated as follows: 



jets 
'-IP 



^ n 



kinematic region defined in Tab. |1] as a function 
of many variables ||2^. Here, only a selection of 



with Mjj the invariant mass of the dijet system. 
Differential cross sections are measured in the 



the available distributions is presented. 

Figure ^ shows that applying the results of 
Regge-QCD fits, as obtained from the diffractive 
structure function [Q, works very well. The anal- 
ysis of the inclusive diffractive cross section yields 
two different fits with similar x^/ndf (cf. inset on 
top of fig. 1^), each with very similar quark den- 
sities but one with a fiat gluon density function 
("HI fit 2") and another with a gluon density 
function that peaks at high values of the frac- 
tional momentum z ("HI fit 3"). For the dijet 
cross section analysis it seems that the fiat gluon 
distribution is slightly preferred over the peaked 
gluon density. Including the contribution of re- 
solved photons further improves the agreement 
with the data, as can be expected, because for 
a large fraction of the dijet event sample the jet 
transverse momentum squared is actually larger 
than the virtuality of the photon (p|. > Q"^). 

The dijet event sample also provides an inde- 
pendent cross-check of the measurement of the 
pomeron intercept ajp(O), as obtained from the 
Regge-QCD fits to the inclusive cross section. Fig- 
ure ^ displays the xjp distribution of which the 
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Fig. 4. The differential dijet cross section is dis- 
played as a function of zjp and compared to differ- 
ent soft colour neutralisations models. The full and 
dashed-dotted (blue) lines represent the old and new 
(with general area law) soft colour interaction model; 
the dashed (red) line represents the semi-classical 
model. 



slope is directly related to ajp(O). Due to the dijet 
requirement, the kinematic constraints are rather 
strict and the distribution is seen to rise instead of 
fall with increasing xjp. The distribution is, how- 
ever, still very sensitive to a change in the value 
of the pomeron intercept, as can be seen from 
the different Monte Carlo curves corresponding 
to a range in ajp(O) between 1.08 and 1.4. Fig- 
ure ^ also shows the contribution from secondary 
reggeon exchanges which is small and only visible 
at large X]p. Note that the beta distribution ex- 
tends to much lower values than in the analysis of 
the inclusive cross section, where /3-values down 
to 0.04 are reached. 

Figure ^ compares the same experimental data 
to several soft colour neutralisation models. The 



original soft colour interaction model |T5[ and the 
semi-classical model [|TU| roughly agree and show 
the same trend as in the experimental data, al- 
though their normalisation is too low by a fac- 
tor of two. The new soft colour interaction model 
(implementing the generalised area law) has 



Fig. 5. The differential dijet cross section is dis- 
played as a function of zjp and compared to the 
saturation and two-gluon models. The requirement 
X]p < 0.01 has been applied in additional to the 
kinematic cuts of Tab. |l|. The full, dashed and dot- 
ted (blue) lines are predictions from the two-gluon 
exchange model for the qqg component with a 
cutoff at 1 GeV'^, the qqg component with a cut- 
off at 0.5 GeV^ and the qq component, respectively. 
The dashed-dotted (red) line represents the satura- 
tion model. 



good normalisation, but the shape of the zjp dis- 
tribution does not agree with the data. 

Finally, Fig. ^ compares the Z]p distribution 
obtained from data with the saturation model 



and two-gluon exchange model Here an ad- 
ditional cut on xjp has been applied to avoid the 
valence quark region and to exclude sub-leading 
exchanges. The saturation model is somewhat low 
while the two-gluon exchange model roughly de- 
scribes the data. As shown, the qq component 
of the two-gluon model is nearly negligible and 
the cross section is by far dominated by qqg par- 
tonic states. The two-gluon model requires a cut- 
off at low transverse momenta of the gluon in qqg 
states. Lowering this cutoff to 0.5 significantly in- 
creases the cross section and affects mostly the 
low pt pomeron remnant. The model is not able 
to describe data at xjp > 0.01. 
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HI 


ZETIS 


z <^ ^/ ^ luu ^e V 




05 < ?/ < 7 


02 < ?/ < 7 


Pt{D*) > 2 GeV 


1.5 < prP*) < 8 GeV 


\7]{D*)\ < 1.5 


< 1-5 


xp < 0.04 


xp < 0.016 


My < 1.6 GeV 


/3 < 0.8 


|t| < 1 GeV^ 





Table 2. Kinematic domains used by the HI and 
ZEUS collaborations in the analysis of open charm 
production in diffractive deep inelastic scattering. 

DifFractive open charm production 

The production of charm quarks in diffractive 
processes provides another test of the underlying 
dynamics of diffraction models |2^. Here, it is the 
large charm quark mass that makes QCD calcula- 
tions possible. Moreover, open charm production 
is free of uncertainties relating to bound states, 
like in exclusive heavy vector meson production 
(e.g. ep eJ/ipp)- Just as for dijets, charm pro- 
duction is highly sensitive to the role of gluons in 
the diffractive exchange. 

Open charm production is studied through the 
production of D*^ mesons. Two different decay 
channels have been used: 



TT 



K-71+; 



(9) 
(10) 



with branching ratios of resp. 2.63% and 5.19%. 
Here and in the following the charge conjugate 
channels are always implied. 

The ZEUS Collaboration obtained preliminary 
results using an integrated luminosity of 44 pb~ , 
yielding 85 D* candidates. The invariant mass 
spectrum is shown in Fig. The cross section is 
measured in the kinematic domain listed in Tab. 
and amounts to: 



ZEUS{prel.)/ 



{ep eD*^Xp) = 281 ± 41+^^ pb 



Figure ^ shows the differential cross as a func- 
tion of various observables. In general, the par- 
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Fig. 6 

the form a{AM 
D*^ candidates. 



slow' 



M{K Tr~^). The curve is a fit of 
- M^+)^ + Gaussian, yielding 85 



tonic pomeron and two-gluon models give a rea- 
sonable description of the data. From the P spec- 
trum however, it is clear that one needs to com- 
bine the qq and qqg components of the two-gluon 
model in order to describe the full distribution. 
The partonic pomeron model seems to yield a (3 
distribution that is too steep. 

HI also obtained preliminary results on D*^ 
production based on an integrated luminosity of 
21 pb^^ p6|. Due to differences in detector accep- 
tances, the kinematic domain (listed in Tab. 0) of 
the cross section measurement is, however, quite 
different from the one used by ZEUS, making 
direct comparison between the two experiments 
difficult. The cross section as measured by HI 
amounts to: 



a 



Hl{prel.) 



{ep eD*^XY) 



154 ± 40 ± 35 pb 
(12) 



Although the differential cross sections obtained 
by HI and the models agree in shape, here, the 
normalisation of the overall cross section is very 
different. 

Because of the different conclusion reached by 
HI and ZEUS, the need for a comparison be- 



8 



Pierre Van Mechelen: Final States in Diffraction at HERA 



1 



ZEUS 1995 - 1997 Preliminary 

:^ 3 



10 



10 



C^'(GeV') 



"L 



50 100 150 200 250 

W(GeV) 



s 

o 

■§200- 



q> 

§ 

I 



10- 



10 



"^2000 



0.005 0.01 



0.015 

y 

pom 



1500 
1000 
500 




0.2 0.4 0.6 0.8 



.300 



200 



100 







• ZEUS 1995-1997 Data 






— ACTWFITD 








— RIDIqq 








-- RIDIqq g 






























1 


1 


























I 


( 









Pj{D) (Gevj 



^(D) 



Fig. 7. Differential 
cross sections for the 
reaction ep — > eD*^Xp 
are plotted as a function 
of various observables. 
The D*^ pseudo- 
rapidity r]^ and 



transverse 



momentum 
are defined rela- 
tive to the 7* axis in 
the rest frame of X. 
The full (black) and 
dashed (blue) lines are 
predictions from the 
two-gluon exchange 
model and represent 
the qq and qqg com- 
ponent, respectively. 
The dotted (red) curve 
is the prediction of 
the partonic pomeron 
model where the flat 
gluon distribution ("HI 
fit 2") is used as input. 



tween experimental results is evident. The RAP- 
GAP Monte Carlo has been used to make pre- 
dictions for the different kinematic domains used 
by the experiments. The result listed in Tab. ^ 
are obtained using the partonic pomeron model 
with "HI fit 2", a charm quark mass of rric = 
1.35 GeV, Aqcd = 0.25, 5 flavours and the evo- 



lution scale set to /i^ = Q"^ + Px + 4:ml. The 
contribution of sub-leading reggeon exchange is 
added separately. The RAP GAP predictions are 
in agreement with the ZEUS results, but are a 
factor three times larger than the HI result. The 
question remains whether this can explained by 
the inability of the model to bridge the gap be- 
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in 


ZEUS 


experiment 
RG P only 

RG p + m. 


154 ± 40 ± 35 pb 
466 ± 12 pb 
504 ± 15 pb 


281 ± 4ll^^ pb 
309 ± 9 pb 
320 ± 12 pb 



Table 3. Comparison of the measured and predicted 
cross sections in the kinematic domains as used by HI 
and ZEUS. The prediction are based on the RAP- 
GAP Monte Carlo model with parameters as ex- 
plained in the text. 
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tween kinematic domains. It is clear that new and 
more precise experimental data are necessary to 
resolve this issue. 

5 Summary 

Recent preliminary results, obtained by the HI 
and ZEUS collaborations, on high Et dijet pro- 
duction and open charm production in diffractive 
deep inelastic scattering have been presented and 
compared to several theoretical models. 

The partonic pomeron model, with a flat gluon 
density function ("HI fit 2"), produces a very 
good description of all experimental data. The 
dijet analysis in particular thus confirms the pic- 
ture of a hard pomeron (q;jp(0) = 1.2) dominated 
by gluons carrying large fractional momenta. 

The soft colour interaction, semi-classical and 
saturation models are either wrong in normalisa- 
tion or in shape. As some of these models might 
be improved by including next-to-leading order 
contributions, it is presently not yet possible to 
make any final judgements on their validity. 

Full perturbativc QCD calculations, based on 
the exchange of two gluons in a colour-singlet con- 
figuration, are able to describe the differential di- 
jet cross sections, but are limited to a kinematic 
domain where xjp < 0.01. 

HI and ZEUS results on open charm produc- 
tion in diffractive deep inelastic scattering suffer 
from low statistics and are, at the moment, diffi- 
cult to reconcile. Until new and more precise data 
become available, it is not possible to draw any 
firm conclusions based on these results. 
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